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Abstract 

A method is proposed for the determination of the position and inclination 
angles of the plane of a spiral galaxy based on the assumption that every spiral 
arm is a monotonia function of the radius versus azimuthal angle. This method 
may yield more accurate results than the more commonly employed isophote 
method, which is fraught with various drawbacks. The use of the new method 
is illustrated by applying it to a sample of 43 objects, and the results agree well 
with data from other sources. 



1. INTRODUCTION 

The orientation of a galactic disk with respect to the line-of-sight is one of the most 
important parameters of a spiral galaxy required for studies of its morphology and dynamics. 
This orientation is determined by the position angle PA and inclination angle i of the plane 
of the disk relative to the plane of the sky. Several methods for determining these angles 
are known, some making use of kinematic data and other based on photometry. The latter 
category includes the most widely used method for determining the inclination angle i of a 
galactic plane relative to the plane of the sky — the so-called isophote method. This method 
estimates i from the degree of elongation of the observed elliptical galaxy isophotes, with or 
without allowance for the finite thickness of the disk. However, the results obtained with 
this method have large uncertainties if the brightness distribution in the galaxy deviates 
strongly from axial symmetry. In fact, this is the case in many observed galaxies, due to 
the presence of bright spirals, a bar, and/or dust lanes; for isophotes in regions of the disk 
at large galactocentric distances, warping of the disk can also contribute to deviations from 
axial symmetry. The method that we propose here can yield more accurate results in such 
cases. 
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All the methods employed assume that the galactic disk can be described using some 
model. Our method for determining the angles PA and i is based on a simple and very 
general assumption about the shape of every spiral arm. We assume the spiral arm to be 
monotonic; i.e., the radius vector R of the spiral in the galactic plane increases monotonically 
with 0: ^ > 0. Section [2] describes the main features of the method and justifies the adopted 
assumptions. Section [3] describes analytically the capabilities of the method as applied to 
logarithmic spirals, which are simplest in terms of the mathematics involved. Section [4] 
compares the method with other methods and discusses its advantages and disadvantages, 
while Section [5] presents the results of applying the new method to a sample of 43 galaxes 
from the Ohio State University Bright Spiral Galaxy Survey. We present the derivation of 
all formulae used in the Appendix. 



2. IDEA BEHIND THE METHOD 

The method we describe here is based on the intuitively clear assumption that all spiral 
arms except for ring-shaped structures are monotonic. In other words, the radius increases 
monotonically in the course of moving along a spiral from the center of the galaxy toward 
its periphery. Note that this should be true not only for ideal, two-armed (or many-armed) 
spirals, but also for branching spirals. Spirals can be traced using young stars (in B images), 
gas (Ha and Nil images), or even dust lanes (5, F, and R images). In this case, the old 
stellar population is of little use due to the relatively large width of the corresponding spiral 
arms. 

The method proposed is not trivial. As is shown in Fig. [1], when a spiral ia a mono- 
tonic function in the plane of the galaxy — i.e., the radius vector R of the spiral increases 
monotonically with the angle 6 — the same spiral may obey a non-monotonic function in the 
plane of the sky. Suppose that we choose the center of the galaxy as our coordinate origin, 
place several points with coordinates and ipi (the radius and azimuth in the plane of the 
sky) onto the spiral arm, and calculate for various PA values ranging from 0° to 180° and i 
values ranging from 0° to 90° corresponding coordinates -Rj(rj, PA, i) and Oiiri, ipi, PA, i) 
in the deprojected plane; the values of PA and i for which the spiral is monotonic (i.e., 
9i{Ri) is monotonic) then give the ranges of possible inclination and position angles for the 
given spiral. 
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Fig. 1. — Schematic demonstrating how the monotonic behavior of the functions R{0) and 
r{ip) differs in the planes of the galaxy and of the sky. Whereas R{6) is monotonic in the 
plane of the galaxy ^ > , r{ip) may be non-monotonic; i.e., ^ may change its sign. This 
is especially obvious in the case of small pitch angles and large inclination angles i. Note 
that we use the polar coordinates i?, 6 and r, ip in the planes of the galaxy and of the sky, 
respectively. 



3. TESTING THE METHOD 

We will first analyze the capabilities of the method described, which we will refer to as 
the monotonic spiral arm (MSA) method. As we pointed out above, this method can be 
applied to any spiral. However, to simplify the calculations, we analyze theoretically how it 
operates with a logarithmic spiral, 6 = p\nR, where p determines the pitch angle, defined 
as the angle between the tangent line to the spiral curve and a perpendicular to the radius 
vector at the given point. It is easy to show that, given the equation of the curve, the the 
tangent of the pitch angle is given by = ^. It is obvious that the logarithmic spiral 
0{R) is a monotonic function. 

Let us first note that a logarothmic spiral with pitch parameter p inclined at an angle 

i smaller than i^r = arctan "^^^^^^^^^ iKl appear monotonic even in the plane of the sky (see 
the Appendix 18.11 for the derivation of this statement). Thus, for a face-on spiral (i.e., a 
spiral lying in the plane of the sky) the method described yields the range of inclination 
angles between to • It can easily be shown that, for a spiral inclined by an arbitrary 
angle i, our method yields PAmsa and ^msa values in the intervals [see the Appendix for the 
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Fig. 2. — Widths of the intervals of PA (left) and i (right) determined using the MSA 
method for logarithmic spirals as functions of the true inclination angle i. The right plot 
shows two curves for each p: Ai~ (the lower curve) and Az"*" (the upper curve). 



derivation of (fTSll and (l25ll]: 
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It is obvious that the true i and PA values are located within these intervals. Figure [2] shows 
the widths of the resulting intervals of «msa and PAmsa as functions of i (it is clear that the 
PA does not affect the result, since a change in PA reduces to a simple turn of the entire 
system around the line-of-sight). 

To apply the method described we wrote a code, which allows to select the center and 
indicate points onto spiral arms. The output is a drawing of angles found on PA-i diagram. 
The program is available upon request (poltorak@alsenet.com). 

We first applied our method to a simulated, but realistic galaxy image with logarithmic 
spiral arms. We adopted for the dependence of the brightness of points on the galactic 
coordinates the relation 



2 

cos(9-plni;) ' 
sR 
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where A, Rq, Ri,a, s, and p are constants of the modeU. We then inclined the simulated 
galaxy using different angles PA and i. 

Applying our method to these images demonstrated that: 

- the method yields a domain of possible orientation angles for the simulated galaxy that 
contains the true PA and i values; 

- Ai and APA are in excellent agreement with the results of analytic computations 
[formulae [U and [2] ; 

- the MSA method is not very sensitive to the number of points on the arm (we rec- 
ommend separating the points by at least 10 times the pixel size in order to reduce 
discretization errors); 

- to determine the angles PA and it is sufficient to follow a half-turn of the spiral with 
six uniformly spaced points; 

- analyzing different parts of the same or several arms may improve the accuracy and 
indicate the presence of possible disk warps at the periphery (manifest as a change in 
PA and i); 

- the uncertainty introduced by inaccurate choice of the center is usually significantly 
smaller than the derived widths of the angle intervals; 

- the MSA methods works well for galaxies with large inclination angles and tightly 
wound spirals (i.e., with large p), but is essentially useless for almost face-on galaxies. 



4. COMPARISON WITH OTHER METHODS 

Having described the MSA method, we now compare it with other methods for deter- 
mining the orientation angles of galaxies. 



^ Here, we put meanings A — 0.7, Rq = Ri = 100, a = 30, s — 0.001, and p = 2 7. The image was 
500 X 500 pixels in size and the galactic center was located at the center of the image. We included the 
term to keep the width of the spiral arms approximately constant, while the term sR slightly broadens 
the arms toward the periphery. 
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Photometric Methods 
Fitting elHpses to galaxy isophotes. 

This method assumes that the galactic disk is axially symmetric, and therefore works 
well for early-type galaxies or galaxies with low-contrast spir als. Since thi s method is 
sensitive to the presence of a bar or spiral structure (see, e.g., Istockl kmd )). it can be 



applied only to the outer parts of galaxies, where many galaxies have warp and hence 
where this method is not very accurate. 

Fitting logarithmic or other forms of spirals to the observed spiral arms. 

This method can be applied to galaxies with well-developed spirals, i.e., late-type 
galaxies. However, logarithmic spirals do not describe all spirals well, especially, since 
the spiral arms in many galaxies branch, and often become open at the disk periphery. 
Like our method, this technique implicitly assumes that the spirals are monoto nic, but 



i mpos es stricter limitations. The first version of this method was proposed by iDanver 



(ll942|). 



Maximizing the axisymmetric component of the Fourier expansion. 

This method is based on maximizing the zeroth or minimizing the second coefficient of 
the azimuthal Fourier expansion of the deprojected image introduced by the inclination; 
i.e., this method likewise assumes that the disk is axisymmetric. It also assumes that 
the second harmonic is only due to the disk inclination, which, naturally, is n o t true 



when the galaxy has a two-armed spiral pattern. However, as iBarbera et al.l (120041 ) 
showed for the BAG method described in their paper, this method is fairly reliable, 
even in the presence of a prominent spiral pattern. 



Kinematic Methods 



Approximation of the velocity field using a model with purely circular ro- 
tation. 

The only widely used kinematic method for determining the orientation of a galaxy 
consists in fitting the two-dimensional observed velocity field w ith the velocity field o f 
a model with pur e ly cir cular motions. Howeyer, Stock (jlQSsI ). Fridman et al. ( 1997 ). 
Lyakhovich et al. ( 1997 ). and Fridman et al. ( 20051 ) showed that the presence of a bar 
or spirals introduces systematic deviations into the derived PA and i values due to 
presence of non-circular motions. This method also cannot be applied to large galaxy 
samples, since it requires sophisticated measurements with a Fabry-Perot interferom- 
eter. 
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The Proposed Method 

The MSA method can be classified as a photometric method. However, unhke other 
methods, it is based on a visual processing of the image: the observer must choose points 
on the spiral arm. This approach has its disadvantages: an astronomer with an active 
imagination could obtain incorrect results by choosing the trace points to be in the regions 
of low signal-to-noise ratio. The advantages of the method include the fact that it does not 
require the image to be prepared for reduction (i.e., there is no need to mask foreground 
stars, subtract the sky background, etc.), that makes this method easy and fast to apply to 
large samples. As we mentioned above, our method is most accurate when applied to highly 
inclined galaxy disks with well defined tightly wound spiral structure. Note also that, unlike 
other methods (except fitting logarithmic spirals to the observed arms), our technique makes 
use of the shape of the spiral when determining the orientation angles. 



APPLICATION OF THE METHOD TO THE OHIO SAMPLE 



To test our method on a large sample of real ojects, we applied it to the 43 galaxies 
from the Sta te University Bright S piral Galaxy Survey, and compared the results with those 
reported by iGarcia-Gornez et al.l (120041) and given in the HyperLeda catalog. Our results 
agree well with those of iGarcia-Gomez et al.l (l2004l ). obtained using the BAG method. The 
agreement with the HyperLeda data is slightly poorer, but this can be explained by the 
relatively modest quality of the data in this catalog (the data were adopted from various 
sources and are of various qual ities; some of the data were adopted from photometric mea- 
surements in the RC3 catalog (Ide Vaucouleurs et al.l (Il99ll )). etc.). Figure [3] compares our 
PA and i values with PAbag and PAnyperLeda, and isAG and i HyperLeda- Figure [4] shows the 
errors obtained using our method as a function of the inclination angle i. As expected, the 
widths of the resulting domains depend on the inclination angle, but are smaller than those 
calculated in Section [3l This may be due to the small disk warps involved: we used various 
parts of all the spiral arms, and our method yielded the intersection of the domains for the 
individual spiral parts. This conclusion is supported by the fact that small (half-turn) parts 
of the spirals yield angle domains that are wider and somewhat different from those derived 
using the entire spiral. 

As an example. Fig. [5] shows an image of the galaxy NGC 1300 with the tracing points 
superimposed, and Fig. [6] shows the results of applying the method to this galaxy. A list 
of all galaxies whose orientation angles relative to the line-of-sight have been determined is 
presented in Table [H 
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Fig. 3. — Co mparison of our calcu l ated values of PA (left plots) and i (right plots) with 
those results of lGarcia-Gomez et alJ (120041 ) using their BAG method (upper plots) and those 
data from the HyperLeda catalog (lower plots). A 5° corridor is shown on both sides in all 
plots. 



6. FURTHER DEVELOPMENT OF THE IDEA 

We propose to further develop this approach by adding an expansion of the image into 
a Fourier series in angle. Note that a correct choice of the orientation angles would make 
the second harmonic of the deprojected image run exactly along the spiral arms. This would 
enable us to refine the angles PA and i. It goes without saying that this proposed improved 
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Fig. 4. — Widths of the derived PAmsa (left) and ^msa (right) intervals as functions of 
^MSA (dots). Also shown are the calculated widths for logarithmic spirals with various p 
values. We consider the MSA method to be inapplicable to galaxies with wide allowed angle 
intervals, and it is this that gives rise to the absence of points in the upper and right parts 
of the plots. 



version of the method can work only for galaxies with symmetrical spirals. Automating this 
method is a fairly complicated task, but, in some special cases, the angles required to meet 
the criterion described can be chosen by hand. 



7. ACKNOWLEDGMENTS 

We thank V.L. Afanas'ev, A. A. Boyarchuk, S.N. Dodonov for fruitful discussions and 
A.V. Zasov for critical comments. This work was partially supported by the Russian Foun- 
dation for Basic Research (project code 05-02-17874-a), the Program of Support for Leading 
Scientific Schools of the Russian Federation (grant NSh. 7629.2006.2), and the Basic Re- 
search Program of the Department of Physical Sciences of the Russian Academy of Sciences 
„Extended objects in the Universe". This research has made use of data of the Bright Spiral 
Galaxy Survey prepared by researchers from Ohio State University, which is supported by 
the National Science Foundation (grants AST-9217716 and AST-9617006) and Ohio State 
Universe. 



- 10 - 



8. APPENDIX 
DERIVATION OF FORMULAE 

Here, we derive the formulas for the domains of possible orientation angles in the case 
of a logarithmic spiral. 



8.1. Critical Inclination Angle 

Let us now derive a formula for the critical inclination angle icr, such that, if tilted by 
no more than this angle a logarithmic spiral remains monotonic. We use two coordinate 
systems: polar coordinates r and in the plane of the sky and R and 6 in the plane of the 
galaxy. Let us now rotate the coordinate systems so that their x axes coincide with each 
other and with the position angles (Fig. [T]). 

We now consider a logarithmic spiral in the galactic plane: 

^ = plni? + ^o or i? = exp l^^-^y (3) 

The coordinates r and ip can be written in terms of R and 6: 

= i?2(l - sin^^ sin^i), (4) 
tany9 = tan6' cosi. (5) 

We now derive an equation for the projection of the logarithmic spiral onto the plane of the 
sky. We begin by using (l5|) to write sin^ 6* as a function of the sky-plane coordinates: 

• 2/) tanV 

^ = — 2-n- — ~- ^ 

cos"^ I + tan^ Lp 

We now substitute ^ and ([6]) into ([4]) to obtain after simple trigonometric manipulations 

2 1 + tan^ if ( 1 N /I n\ / X 

r = — 5 — exp -farctanfA tany9) — b'o) , (7) 

1 + A"' tan \p / 

where we introduced the quantity A = — 

^ cos I 

Since r > by definition, the condition that the spiral increase monotonically, > 0, 
needed to determine (or Acr) can be replaced by the condition 

d(r^) dr , , 
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We now introduce the derivative we need to derive: 

d(r^) „ 1 + tanV ( ^ 1 / 1 , 

-27^^T— i tanV ^ + tan(/? — -1 



dip (1 + A2^tan2v9)2 \ A^p V-^ct / -^crP 

X exp ^-(arctan(Acrtanv3) — ^^o)^ • (9) 

This quantity is greater than zero for any (p if the discriminant of the quadratic equation 
in curly braces (an equation in tan ip) is negative: 

1 :\' 4 



p^Xt,-2i2 + p')\l+p' <0. (11) 
And we finally derive the condition for monotony of the projected, inclined logarithmic spiral: 

1 + ^(1 - ^/^Tl) < Xl < 1 + ^(1 + ^/^Tl). (12) 

By definition, 

-^cr = — T— = 1 + tan^ i„ ^ 1, (13) 

cos"' «cr 

which means that the left-hand inequality in (fT2l l is automatically satisfied for all Acr- Since 
0° ^ ^ 90° we obtain 



'2 (l + Vp' + I 
P 

Hence, the inclinations of the galactic plane relative to the plane of the sky for which the ob 



0^tanzcr<-'^ • (14) 

P 



server will see the spiral as a monotonic function lie in the interval from to arctan 



8.2. Range of Inclination Angles 

As for a face-on spiral, in the case of an arbitrarily inclined spiral, our method yields 
a range of possible inclination angles in the vicinity of the true value. To find this domain, 
let us consider a spiral inclined at an angle i to the plane of the sky. We seek the range of 
inclination angles ^msa such that deprojection yields a monotonic spiral. In addition to the 
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coordinate systems used above, we introduce another coordinate system {R, 9) in the plane 
onto which we deproject the galaxy using the angle ^msa^ 

= R\l - sin^ 9 sin^ ^msa), (15) 
tan(y9 = tan6' COS2MSA- (16) 

We obtain, similarly to what we suggested in the preceding clause: 

-2 l + tan^fii f2 



R = —r=exp [ -(arctan(Atan6') — ^^o) , (17) 

l + \'^ta.ir9 \P J 

where we have changed A: A = ^°cosT^ • order for the deprojected spiral to be monotonic, 
it is necessary and sufficient that the derivative '^^^ •* not change sign. Noting that (I7j) and 
f flTll are identical in form, we can immediately write the answer: 



'i + 4(1 - < "-^^^ < Ji + 4(1 + v^). (18) 

cost V P 

This is the final formula for the desired range of inclination angles yielded by our method. 



8.3. Range of Position Angles 

To find the range of position angles, we consider the spiral deprojected using the exact 
inclination angle i and various position angles PAmsa = + APA. 

Here, we introduce another supplementary coordinate system {R, 9) in the plane rotated 
through the angles i and PAmsa- We align the x axis of this system along the assumed 
position axis, i.e., we rotate it by an angle APA relative to the x axis in the plane of the sky 
(which coincides with the direction of the true position axis) about the line-of-sight. The 
introduced coordinates can be related to the coordinates in the sky plane as follows: 

= p2(i_sin2^sin2z), (19) 

cosAPA tan^ + sin APA/cosi , , 

tan 09= = '— . (20) 

-sin APA tan6' + cos APA/cosi 



We similarly derive 



R' = f{9)R' = f{9) exp {^-{g{9) - ^o)) , (21) 
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where 

1 + tan^ 9 cos^ i + sin^ i (tan 9 cos APA + sin APA / cos z) ^ ' 



Aan(v9(6'))\ / 1 tan 6* cos APA + sin APA/ cos 2 ' 

g[9) = arctan = arctan = ; 

\ cosi / \ cos 2 —tan sin APA + cos APA/ cos 2^ 



It can easily be shown that 



it follows that 



= Imm- (22) 

The derivative of f{9) equals 

df(9) 2p(9) sinAPA cosAPA sin^i / „- - ^ r. ^ f ■ 1 

= — i- tan^6' + tan6' tan APA cos«H 

d9 1 + tan^ 9 cos i \ \ cos i 

As a result, we obtain 

am _ d(R^)^^^^ ^ ^,cVg) ^ 41 f 2^.(,-) ^ m\ . (23) 



d9 ' d9 f{9) [p-' ' ' d9 j 

This last quantity is greater than zero if and only if the expression in curly braces is positive: 



tan^ 9 ( + sin^ i sin APA cos APA ) + 
P 

+ tan 9 sin^ i sin^ APA ( cos i 



cosz 



(COS 2 \ 
sin^ i sin APA cos APA J > 0. (24) 

The inequality (l24l l can be satisfied for any 9 value only if the discriminant of the quadratic 
equation (in tan^) is negative: 

sin"^ APA sin^ i ( cos i H : ) +41 sin^ APA cos^ APA sinS-^ <0 

cos 2 J \ 



or 



. . „ , , cos I 
smAPA < 



sin^ z 



2 ( -1^-1). (25) 
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Hence, we find that, when applied to a logarithmic spiral with pitch parameter p inclined 
at an angle i, our method yields the range of position-angle values near the true value from 
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Table 1: The complete list of 43 galaxies whose PA and i have been determined using the 
MSA method. 
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Fig. 5. — Application of the MSA method to the Southern arm of NGC 1300: the B image 
with arm-tracing crosses (one is white to make it easier to see) . 
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Fig. 6. — i-PA diagrams obtained by applying the MSA method to NGC 1300. Left: 
result obtained by applying the method to the Southern arm only. Right: result obtained 
by applying the method to both the Northern and Southern arms. The black domain shows 
the intersection of the gray domains obtained for the individual arms. It is clear that both 
arms yield close values for PA and i. The region of intersection makes it possible to reduce 
the intervals of possible values. 




